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In the rotational sandpile model, either the clockwise or the anti-clockwise toppling rule is assigned 
to all the lattice sites. It has all the features of a stochastic sandpile model but belongs to a different 
universality class than the Manna class. A crossover from rotational to Manna universality class 
is studied by constructing a random rotational sandpile model and assigning randomly clockwise 
and anti-clockwise rotational toppling rules to the lattice sites. The steady state and the respective 
critical behaviour of the present model are found to have a strong and continuous dependence on 
the fraction of the lattice sites having the anti-clockwise (or clockwise) rotational toppling rule. As 
the anti-clockwise and clockwise toppling rules exist in equal proportions, it is found that the model 
reproduces critical behaviour of the Manna model. It is then further evidence of the existence of 
the Manna class, in contradiction with some recent observations of the non-existence of the Manna 
class. 

PACS numbers: 89.75.-k,05.65.+b,64.60.av,68.35.Ct 


I. INTRODUCTION 

A sandpile is a prototypical model to study self- 
organized criticality (SOC) [l[, which refers to the in¬ 
trinsic tendency of a wide class of slowly driven sys¬ 
tems to evolve spontaneously to a non equilibrium steady 
state characterized by long-range spatiotemporal correla¬ 
tion and power-law scaling behaviour. Several crossover 
phenomena from one sandpile universality class to the 
other are reported in the literature on sandpile models. 
For example, a crossover from Bak, Tang and Wiesenfcld 
(BTW) [2j to the stochastic Zhang model was observed 
by O. Biham et al. || by controlling the fraction of en¬ 
ergy distributed to the nearest neighbours in a toppling. 
A crossover from the deterministic Zhang model [J] to 
the stochastic sandpile model (SSM) [5H3] was studied 
by Liibeck l§| by controlling the threshold condition. A 
crossover from the directed sandpile model (DSM) Q 
to the directed percolation (DP) class 0 was observed 
by Tadic and Dhar by introducing a stickiness parame¬ 
ter in the DSM ET The crossover phenomena studied 
in these models are usually from a deterministic to a 
stochastic model. However, the universality class of a 
sandpile model is believed to be determined by the un¬ 
derlying symmetry present in the model |12 :.]. A crossover 
from one sandpile universality class to another then re¬ 
quires a change in the underlying symmetry of a given 
model. It is therefore intriguing to study a crossover 
phenomenon within the stochastic class of models with 
different symmetries in the toppling rule, and to look 
for spontaneous symmetry breaking in the system as one 
of the system parameters is tuned. Two such stochas¬ 
tic sandpile models are the SSM [7j and the rotational 
sandpile model (RSM) jl3j, 0. The SSM is governed 
by externally imposed stochastic toppling rules. On the 


other hand, the RSM is governed by deterministic rota¬ 
tional toppling rules (except the very first toppling) and 
has broken mirror symmetry. Such a model can be useful 
in studying the avalanche dynamics of charged particles 
in the presence of a uniform magnetic field. In the RSM, 
an internal stochasticity appears due to a superposition 
of toppling waves from different directions during time 
evolution. Eventually, that induces all the features of a 
stochastic sandpile model, such as toppling imbalance, 
negative time auto correlation, and existence of finite- 
size scaling (FSS) into the RSM 0-0. The RSM is 
thus a stochastic model, but it belongs to a completely 
different universality class than the Manna class of the 
SSM. The question is whether it would be possible to 
reproduce the critical behaviour of the SSM of Manna 
type in a model such as the RSM, which is stochastic 
due to its internal dynamics. Moreover, there is a long 
standing debate in the study of SOC as absorbing state 
phase transitions (APT) 0 - 0 ] stating that the stochas¬ 
tic universality class or the Manna class is essentially the 
DP universality class. There is continuing evidences in 
favour [2C| and against E] the existence of the so-called 
Manna class. If due to any external condition on the 
RSM, it reproduces the critical behaviour of the Manna 
model, that will be additional independent evidence for 
the existence of the Manna class. 

In this paper, the crossover from one stochastic uni¬ 
versality class to another is studied by constructing a 
random rotational sandpile model (RRSM) and the ex¬ 
istence of the Manna class is discussed in the context 
of random mixing of two conflicting rotational toppling 
rules in the model. 


II. THE MODEL 


The RRSM is defined on a two-dimensional (2D) 
* santra@utg.ac.in square lattice of size L x L. Initially, all lattice sites 
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FIG. 1. The two toppling rules of RRSM, CTR, and ATR are 
demonstrated on a square lattice in (a) and (b), respectively. 
The active site in grey is at the center and it has received 
the last sand grain from the bottom (say), represented by a 
dashed arrow. The directions of sand flow from the active 
site are represented by solid arrows. In (a), due to CTR, one 
sand grain goes to the right and the other goes up, whereas 
in (b), due to ATR, one sand grain goes to the left and the 
other goes up. 


are assigned with the clockwise toppling rule (CTR). A 
fraction p of lattice sites are then changed to the anti¬ 
clockwise toppling rule (ATR) randomly. The toppling 
rules assigned to the lattice sites remain unchanged dur¬ 
ing the time evolution of the system and hence this can 
be considered as a quenched random configuration of the 
toppling rules. The RSM [l3[ was defined for the pres¬ 
ence of only one type of toppling rule either CTR or ATR. 
Since the sandpile dynamics is independent of the sense 
of rotations, the RSMs with either CTR or ATR have the 
same critical behaviour. 



FIG. 2. Plot of h p (n) against n, the number of sand grains 
added, for p — 0.0, 0.1, 0.25 and 0.5 for RRSM on L — 2048. 
In the inset, the average steady state height ( h(p )) is plotted 
against p. The dashed line represents the average steady state 
height of the SSM studied on L = 2048. 

by 

CTR : hi —y hi — 2, hj —> hj + 1, j = di , dj+i, (1) 

where one sand grain goes along di and the other goes in 
a clockwise direction with respect to di. If the active site 
i is with ATR, the sand distribution is given by 

ATR : hi — y hi — 2, hj — y hj -hi, j = di , di— i, (2) 


Each lattice site *, irrespective of the type of toppling 
rule it has, is assigned with a positive integer hi repre¬ 
senting the height (the number of sand grains) of the sand 
column. Initially, all /qs are set to zero. The system is 
driven by adding sand grains, one at a time, to randomly 
chosen lattice sites i. The critical height of the model is 
taken as h c = 2. As the height of a sand column hi be¬ 
comes equal to or greater than the critical height h c , i.e., 
hi > h c , the site becomes active and bursts into a top¬ 
pling activity. On the very first toppling of an active site, 
two sand grains are given away to two randomly selected 
nearest neighbours out of the four nearest neighbours on 
a square lattice. As soon as a site j receives a sand grain, 
the direction dj from which the grain was received is as¬ 
signed to it besides incrementing the height of the sand 
column hj by one unit. The value of dj can change from 
1 to 4, as there are four possible directions on a square 
lattice. The directions from an active site i are defined as 
di = 1 for left, di = 2 for up, di = 3 for right and di = 4 
for down. As the avalanche propagates, the direction di 
and height hi are updated upon receiving a sand grain, 
and only the information regarding the direction from 
which the last sand grain was received is kept. The next 
active sites with hi > h c in the avalanche will topple 
following a deterministic rotational toppling rule. The 
toppling rules for an active site i that has received the 
last sand grain from a direction di are given below. If the 
active site i is with CTR, the sand distribution is given 


where one sand grain goes along di and the other goes 
in an anti-clockwise direction with respect to di. If the 
index j becomes 5, it is taken to be 1; if it becomes 0, it 
is taken to be 4. The CTR and ATR are demonstrated in 
Fig. [T| on a square lattice. The avalanche stops if there 
is no active site present and the system becomes under 
critical. The next sand grain is then added. As RSM, the 
RRSM is non-Abelian [23] and it has no toppling balance 


In the following, the results of the RRSM are compared 
with those of the original RSM [l3[ and the SSM 0 • The 
SSM considered here is a modified version of the Manna 
model [S, Q known as the Dhar Abelian model. The 
toppling rule in this SSM is that two sand grains of an 
active site are given to two randomly selected nearest- 
neighbour sites out of four possible nearest neighbours 
on a square lattice and the height of the active site is 
reduced by 2. The remaining sand grains remain at the 
present site. 


III. STEADY STATE 

The steady state of a sandpile model corresponds to 
constant currents of sand influx and outflux. Conse¬ 
quently, the average height of the sand columns remains 
constant over time. For a given value of p in the RRSM, 
the mean height h p (n) of the sand columns is expected 
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to be constant over the number of sand grains n added 
(equivalent to time) to the system. It can be defined as 


h p (n) 


1 

L2 


L 2 

^ ~2 h i(p,n ) 

i=i 


(3) 


for a system size L. For different values of p, h p {n) for 
L = 2048 is plotted against n, the number of sand grains 
added in Fig. [2j After the addition of a sufficiently large 
number of sand grains, the system reaches a steady state 
corresponding to a given value of p. In order to study 
the effect of p on the steady state height, the saturated 
average height (h(p)) of the sand columns in the steady 
state for a given value of p is estimated taking the average 
over the last 10 5 sand grains on every 64 random config¬ 
urations of toppling rules. Note that, no configurational 
average is required for p = 0 and 1. In the inset of Fig. [2j 
( h(p )) is plotted against p for L = 2048. It can be seen 
that the saturated average height of the sand columns 
at the steady state decreases as p is varied from 0 (or 
1) to 0.5, and it attains a minimum value at p = 0.5. 
The values of ( h(p )) are found to be symmetric about 
p = 0.5, as expected. The average heights at p = 0 and 1 
are found to be that of the RSM [l3j, whereas at p = 0.5 
the average height is that of the SSM. The time-averaged 
steady state-height for the SSM for L = 2048 is measured 
independently, and it is found to be 0.7162 ± 0.0001. It 
is represented by a dashed line in the inset of Fig. [2] It 
can be noted here that the measured value of the aver¬ 
age sand column height for the SSM is in good agreement 
with that of the driven dissipative sandpile in the con¬ 
text of the precursor to a catastrophe study Q as well 
as the critical point of the APT of a fixed energy sandpile 
model 0 . Thus, the steady-state heights corresponding 
to different values of p are not only different from each 
other but also very different from that of the SSM. 


IV. RESULTS AND DISCUSSION 


The critical properties of RRSM are characterized 
studying the properties associated with avalanches in the 
steady states at different values of p and the system size 
L on 2D square lattices. The value of p is varied from 
0 to 1, and the system size L is varied from 128 to 2048 
in multiples of 2 for every value of p. For the sake of 
comparison, data for the SSM are also generated for the 
same lattice sizes. The information of an avalanche is 
kept by storing the number of toppling of all the lattice 
sites in an array Sx jP [i], i = 1, • • • , L 2 for given L and p. 
All avalanche properties of interest, such as the two point 
toppling height correlation function, the toppling surface 
width, avalanche size, etc., will be derived from Sl )P [i\. 



FIG. 3. (Colour online) Plot of Sl, p P{Sl, p ) versus Sl, p for 
RRSM at different values of p (0.0, 0.1, 0.5) for L = 2048. For 
the sake of comparison, the same distribution for the SSM is 
shown by a solid line. 



(a) RRSM p=0.50 



(b) SSM 


FIG. 4. (Colour online) (a) Toppling surfaces of RRSM at 
p = 0.5 generated on a square lattice of size L — 128. (b) Top¬ 
pling surface of the SSM generated on the same system size. 
The highest toppling number corresponds to a light brown 
colour and the lower toppling numbers are represented by 
darker and darker brown colours. Below each surface, their 
2D projections are shown. These projections represent the 
avalanche clusters in two dimensions. 


A. Distribution of Sl iP and avalanche morphology 

The probability distribution of Sl, p is defined as 
P{Sl,p) = ns/A where ns is the number of sand columns 
that toppled S times and A is the number of distinct 
sand columns (or lattice sites) toppled. A FSS study of 
P{Sl,p) reported in [25j suggests that RSM and SSM 
follow FSS, but BTW does not. In Fig. [31 distributions 
Sl, p P(Sl, p ) of RRSM are plotted against Sp, P for sev¬ 
eral values of p for a fixed L = 2048 and compared with 
that of the SSM. The distribution S'l, p F > (5'l,p) is found 
to be of Poisson type as expected. However its height, 
width, and the peak position vary strongly with p on a 
given lattice. It can also be noted that the distribution of 
the SSM is not identical with that of RRSM at p = 0.5. 
This implies that the internal structure of an avalanche 
at different values of p is different, but also it is different 
from that of the SSM in comparison to that of RRSM 
at p = 0.5. The avalanche morphologies of two typical 
avalanches generated on a 128 x 128 square lattice for 
RRSM at 0.5 and for that of the SSM are presented in 
Fig. SI The values of the toppling number at different 
lattice sites of an avalanche define a surface in three di- 
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mensions called toppling surface 0 ], Thus, the height 
of the toppling surface at the ith lattice site is then given 
by Sl,p[i}- The toppling surface of RRSM at p = 0.5 
is found to be fluctuating all over the lattice as that of 
the SSM with different maximum heights. For SSM, it 
is approximately 31, whereas that for RRSM at p = 0.5 
is approximately 12, similar to the observation of their 
steady-state heights. The projection of the toppling sur¬ 
faces in two dimensions is shown below the respective 
toppling surfaces. The 2D view of a toppling surface is 
known as an avalanche cluster. It can be seen that the 
avalanche cluster of RRSM at p = 0.5 exhibits random 
mixing of colours representing different toppling numbers 
as that of an avalanche cluster of the SSM. Note that 
both are very different from that of the RSM in which 
a random superposition of several BTW-type concentric 
zones |0 of lower and lower toppling numbers around 
different maximal toppling zones is observed [l.3| . Since 
the avalanche morphologies of RRSM at p = 0.5 and the 
SSM are found to be similar, it is expected that both 
models have the same critical behaviour, though the dis¬ 
tributions Sl, p P{Sl,p ) of them are different. 


B. Properties of avalanche size 


One of the macroscopic avalanche properties is the to¬ 
tal number of toppling in an avalanche, called the top¬ 
pling size s. Knowing the values of S7. p ['i] at every lattice 
site, the toppling size s of an avalanche can be obtained 
as 


L 2 

s(L,p) = J2 Sl ^ ( 4 ) 

i—1 

for given L and p. At the steady state, 5 x 10 5 avalanches 
are generated on every quenched random configuration of 
the toppling rules for a given value of p. For every value 
of p , 64 configurations of quenched toppling rules are con¬ 
sidered. Therefore, ensemble averaging is performed over 
Ntot = 32 x 10 6 avalanches for given values of L and p. 
The probability to have an avalanche of toppling size s 
is given by the Pl, p {s) = N s (L,p)/N tot where N s (L,p) 
is the number of avalanches of toppling size s for given L 
and p out of total number of avalanches N to t generated. 
For RSM (corresponding to p = 0 and 1 of RRSM), it is 
already known that the distribution of s follows a power 
law scaling with a well defined exponent r s and obey FSS 
0,111. The probability distributions Pl, p (s) for the top¬ 
pling size s for several values of p (other than 0 and 1) for 
a fixed large system sizes L = 2048 and for several values 
of L for a fixed value of p show that not only the cutoffs 
but also the slopes of the distributions depend on p for a 
fixed L whereas for a fixed p only the cutoffs of the distri¬ 
butions depend on L keeping the slopes unchanged. For 
a given p , a probability distribution function for toppling 
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FIG. 5. (Colour online) (a) Plot of a s (p, q) against q for differ¬ 
ent values of p. (b) The variations of r s (p) and D a (p) against 
p are shown. 


size s is then proposed as 


Pl, p (*) = s~ Ts(p) fp 


IL d s(p)} ’ 


(5) 


where f p is the p dependent scaling function and D s (p) is 
the capacity dimension of the toppling size s correspond¬ 
ing to given p. To have estimates of the exponents r s (p) 
and D s (p) defined in Eq. ([5]l. the concept of moment 
analysis 23. 0 for the avalanche size s has been em¬ 
ployed. For a given p , the gth moment of s as function 
of L is obtained as 


pOO 

(s q (L,p)) = / s q P L , p (s)ds ~ L 
Jo 


( p,q ) 


( 6 ) 


where the moment scaling exponent cr s ( p , q ) would be 
cr s(p, q) = D s (p)[q+1-T s (p)} for large q as q > t s {jp) - 1. 
For each value of p, a sequence of values of <r s ( p , q) as 
a function of q is determined by estimating the slope of 
the plots of \og(s q (L,p)) versus log(L) for 200 equidistant 
values of q between 0 and 2. cr s (p,q) is plotted against 
q for different values of p in Fig. [51(a) . To obtain the 
values of t „ (p) and D s (p ), the direct method developed 
by Liibeck 27} is employed in which a straight line is fit¬ 
ted through the data points for q > T s {p) — 1. From the 
straight line fitting, the x intercept provides r s (p) — 1 and 
the y-intercept provides _D s (p)[l — r s (p)]. Straight lines 
are fitted through the data points for different values of p 
in the range of q between 1 and 2, and the x and y inter¬ 
cepts are noted. The estimated values of the exponents 
T s (p) and D s (p) are then presented in Fig. [5[b) against 
p. There is a continuous change in the values of r s (p ) 
and D s {p) as p changes from 0 to 1. This indicates a 
continuous crossover of the critical behaviour of the sys¬ 
tem through a series of universality classes of RRSM at 
different values of p. The exponents T s (p) = 1.234(13) 
and D s (p) = 2.82(2) at p = 0 and 1 correspond to that 
of the RSM 0] as expected. However, the values of the 
critical exponents T s {p) = 1.286(9) and D s (p) = 2.71(2) 
at p = 0.5 are found very close to that of the SSM [271 — 
[29j . Although both of the exponents are varying continu¬ 
ously with p , due to the diffusive behaviour of RRSM, the 
scaling relation D s ( p ) [2 — t s (p)] = 2 is found to be valid 
within error bars for all values of p. It can be noted here 
that the continuous crossover in RRSM from one stochas¬ 
tic to another stochastic universality class through a se- 
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FIG. 6. (Colour online) Plot of P L , p (s)L Da(p)ra(p) against 
s/L d “with different values of L is shown in (a) for p = 0.25 
and in (b) for p = 0.50. 


ries of stochastic classes is very different from the ob¬ 
served crossover phenomena from one deterministic to a 
stochastic universality class 3Cj] • 

To verify the form of the scaling function f p given 
in Eq. J5]), the scaled distributions Pl^ p {s)L D s ^ Ts ^ p ' > 
is plotted against the scaled variable s/L Da ^ for four 
different system size L for p = 0.25 in Fig^a) and for 
p = 0.50 in Fig. Ob) in double logarithmic scales. It can 
be seen for both the cases that data for different values 
L are collapsed onto a single curve, ie., the scaling func¬ 
tion. Hence, the proposed scaling function form in Eq. 
© for RRSM for any values of p and L is a correct scaling 
function form. The analysis not only provides estimates 
of the values of the exponents, but it also confirms that 
the model exhibits FSS. The steady-state event distribu¬ 
tion of this slowly driven dynamical system is then found 
to obey power-law scaling at different values of p; RRSM 
then exhibits SOC for any value of p with different sets 
of critical exponents. Moreover, at p = 0.5 of RRSM, the 
appearance of SSM confirms the existence of the Manna 
class. 


C. Properties of toppling surfaces 


The toppling surfaces are obtained for large spanning 
avalanches only. The spanning avalanches are those that 
are touching the opposite sides of the given lattice. For 
a given system size L, a total of N span = 1024 spanning 
avalanches are taken over 64 initial random configura¬ 
tions for each values of p. To study the scaling behaviour 
of a two-point height-height correlation function, the cor¬ 
relation between the toppling numbers of two lattice sites 
separated by a certain distance has to be determined. 
Since the correlation function has to be calculated as a 
function of a continuous variable r, (the distance between 
any two lattice sites), the toppling number of a site is rep¬ 
resented as Sl, p (x) : where x is the position vector of the 
lattice site with respect to the origin of a 2D coordinate 
system instead of using a discrete sequence of toppling 
numbers stored in Sx iP [i]. The square of the difference of 
toppling numbers 8Sl, p (t) at two lattice sites separated 


FIG. 7. (a) Probability distribution of <5Sx, P (r) for RRSM 
with p = 0.10 (O), 0.20 (□), 0.25 (O), 0.35 (A), 0.50 (<) 
with L — 2048 and r = 512. (b) Distribution of 5Sz,, p (r) for 
different L as L = 256 (O), 512 (□), 1024 (O), 2048 (A) for 
fixed values of p — 0.25 and r = 128. 


by a distance r is given by 

SS LiP (r) = \S L ,p(x + r) - S LiP (x)\ 2 , (7) 

where Sl,p(x+t) is the toppling number at x+r for given 
L and p. The probability P[8Sl, p {t )] of a particular value 
of SSL, P (r) occurring for a particular r for a given value 
of L, p is defined as 

P[SS L , p (r)} = (8) 


where ngs L p (r) is the number of pairs of sand columns 
having the desired value of SSL, P {r) and N r is the total 
number of pairs separated by a distance r for all the 
surfaces. To determine Ay, for each surface 100 centers 
are randomly selected. From each center, all possible 
sand columns at a distance r are counted and then added 
for N span = 1024 surfaces. The probability distribution 
P[SSL, P (r)] is then estimated for several values of L , p, 
and r. To guess the form of the distribution function 
-P[&Sx,p(j")]> once is plotted against p for a fixed system 
size L = 2048 and r = 512 and then it is plotted against 
L for a fixed value of p = 0.25 and r = 128 in Fig. 
[7fa) and E^b), respectively. It can be seen that as p 
increases, the cutoffs of the distributions decrease for a 
given system size L. On the other hand, for a given 
p the cutoffs increase as the system size L is increased. 
Hence, following Ref. f3l) the form of the probability 
distribution function P[SS(r,p, L)] is proposed as 


P[SS L , P (r)\ 


r~ 2H (p) [ SS L , p (r)' 

9 Lt r 2ff( P ) ’ 


(9) 


where H(p ) is the p-dependent Hurst exponent, £ is an¬ 
other exponents, and g is the scaling function. 

The correlation between the toppling numbers of two 
sand columns separated by a distance r can be obtained 
by estimating the expectation {5Sl{t,p)). The correla¬ 
tion function for a given L and p, is obtained as 

nOO 

C L , p (r) = / 6S L Ar)P[SSLAr)}d[SS L , p {r)} 

J o 

/»oo 

= r 2H ^ L C ’ / zg(z)dz 

Jo 

~ r 2ff(p) i c (10) 
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FIG. 8. (Colour online) (a) Plot of Il, p {R) against R for p = 
0.0 (O), p = 0.10 (□) andp = 0.50 (A) on a system of size L = 
2048. The straight lines represent the linear least-squares- 
fitted lines through the linear region of the data points, (b) 
Plot of Wl, p against L for the same values of p as in (a) with 
the same symbols. The straight lines are the linear least- 
squares-fitted lines through the data points. In (c) and (d), 
H(p) and x{p) are plotted, respectively, against p. The dashed 
lines in (c) and (d) represent the values of H and \ for the 
SSM. 

where z = SSl, p (t) / (L^r 211 ^) is the scaled variable and 
the value of the integral would be a constant. Notice 
that Cl, p {t) is a system size dependent correlation func¬ 
tion. Such correlation functions also appear in the cases 
of growing interfaces in random media and self-affine frac¬ 
ture surfaces [HJ. In order to determine the Hurst ex¬ 
ponent H(p), and the values of £, integrated correlation 
function Ip p (R ) up to a distance R and the overall sur¬ 
face width Wl, p are estimated. Il, p (R ) and Wl, p are 
obtained as 

Il, p (R)= l C L , p {r)dr ~ R 1+2H ^L^ (11) 

Jo 

and 

Wl, p = C L , p (r)rdr ~ I?^\ (12) 

where x(p) = C/2 + i?(p) is known as the roughness 
exponent. Il, p {R) against R and Wl, p against L are 
plotted in double logarithmic scales in Figs. UK a) and 
lib), respectively for the RRSM at p = 0,0.1 and 0.5 
and for different values of L = 128 — 2048. It can be seen 
that both Il, p (R) an d follow power-law scaling with 
their respective arguments. The slopes are obtained by a 
linear least-squares fit through the data points, and one 
obtains 1 + 2 H(p) from (a) and x(p) from (b). The values 
of the Hurst exponents H (p) and the roughness exponent 
x(p) are plotted against p in Figs. He) andHd), respec¬ 
tively. A few observations are there. First, a continuous 
crossover in the values of the critical exponents has oc¬ 
curred as p changed from 0 (or 1 ) to 0.5. This confirms 



FIG. 9. (Colour online) (a) Plot of Il, p /L 1+2H ^ against L 
for different p. The straight line through data points gives 
the value of £ « 1. (b) Plot of P[5SL, P (r)]r 2H( - p ^ L against 
5Sl, p {v)/L r 2H ^ forp = 0.0, p = 0.1, andp = 0.5 for different 
values of L and r (see the legend). 


the existence of a series of stochastic universality classes 
as observed in the case of an avalanche size distribution 
exponent. Second, not only the values of H(p) « 0.35 
and x(p) ~ 0.82 for the RRSM at p = 0 and 1 is found 
to be the same as those of the RSM [3]) but also the 
values of H(p) ~ 0.21 and xip) ~ 0.70 for the RRSM at 
p = 0.5 are found to be close to those of the SSM Il3 |. 
The dashed lines in Fig. He) andHd) represent the val¬ 
ues of H and x for the SSM. Therefore, the Manna class 
exists in the strong disordered limit of the RRSM. Third, 
comparing the values of \{p) an d H(p) for all values p, it 
is observed that x(p) ~ H{p) + 1/2, which suggests that 
(wl. In that case, the values of x(p) and H(p) obtained 
in the RRSM for different values of p do not satisfy the 
usual Family-Vicsek scaling (no difference in the rough¬ 
ness exponent and the Hurst exponent) [33], rather they 
satisfy an anomalous scaling given by x(p) = C/2 + i?(p) 
with £ = 1 for any value of p. Such a scaling relation 
is already verified for the RSM and the SSM [3 113 in¬ 
dependently. Note that the anomalous scaling resulted 
here due to the system size dependence of the correla¬ 
tion function Cl, p {t) [35| . Finally, the critical exponent 
of toppling surfaces and that of the avalanche size distri¬ 
bution are found to be related as D s (j>) = 2 + x(p) 0 
at each value of p. The results obtained in two different 
methods are then consistent. 

To verify the scaling form of the probability distri¬ 
bution P[SSL, P (r)}, the value of the exponent f must 
be determined. To obtain the numerical value of £, 
Il,p/R 1+2H ^ for R = L are plotted against L for several 
values of p in FigHa). From the slopes of the linear least- 
squares-fitted straight lines, it is found that £ « 1 for all 
values of p, as already predicted by the scaling relations. 
The scaling function form given in Eq. © is now veri- 
fied by plotting a scaled distribution P[5SL,p{r)\r 2H< ' p ^ L 
against a scaled variable SSL, P {r)/Lr 211 ^ for different 
values of p in Fig. H^O- F° r C = A good data collapses 
are observed for different values of L and r at three differ¬ 
ent values of p with the respective values of the Hurst ex¬ 
ponent H(jp) . Thus the proposed scaling form assumed in 
Eq. © is correct. Such distribution functions will then 
be useful to analyze rough surfaces arising in a system 
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with controlled disorder. 


D. Comparison of SSM and RRSM at p = 0.5 

Alhough the macroscopic parameters describing the 
critical states of the RRSM at p = 0.5 and the SSM are 
found to be drastically different, the values of the critical 
exponents are found to be similar. It is then important to 
compare the scaling function forms of both models. The 
probability distributions of toppling size s and that of the 
square difference of toppling numbers SS(r ) are consid¬ 
ered for comparison of their associated scaling functions. 
The mo del-dependent probability distributions for s and 
SS(r) are proposed as 





P s ,m = a m s T3m fm(b m s/L Dam ) (13a) 

and 

Pss,m = c m r~ 2Hm g m [d m 5S(r)/r 2Hm ], (13b) 

where m represents the models RRSM at p = 0.5 and the 
SSM, and a m , b m , c m , and d m are non-universal metric 
factors that contain all non-universal model-dependent 
features such as the lattice structure, the update scheme, 
the type or range of interaction, etc. [36[. The p de¬ 
pendence in both the functions and the L dependence 
in Eq. (113bl) are dropped because they are considered 
either for a fixed value of p or for a fixed value of L. 
Assuming a m , b m and c m , d m all are equal to 1, the 
scaled distributions P s , m s Tsm for several values of L and 
Pss,m , r 2Hrn for several value of r for L = 2048 are plot¬ 
ted in Figs. EH a) and IKJl bi respectively, against their 
respective scaled variables s/L Dsm and 8 S (r) / r 2Hm using 
the same values of r sm , H m and D sm for both models. It 
can be seen that in both cases, the distributions collapse 
onto two different curves corresponding to the RRSM at 
p = 0.5 and the SSM. It seems that scaling functions are 
different for these two models, although they scale inde¬ 
pendently with their respective arguments with the same 
critical exponents. It is then essential to verify whether 
the scaling functions are affected by the non-universal 
metric factors or not. The values of a m are calculated 
from the limiting value of P St mLS Tam as s — > 0, and they 
are found to be a m = 0.34±0.01 for the RRSM at p = 0.5 
and dm = 0.27 ± 0.01 for the SSM. The values of b m 
are calculated from the average toppling numbers s of 
the two models for a fixed L. For L = 2048, the val¬ 
ues are b m = 2 _16 - 92±002 for the RRSM at p = 0.5 and 
b m = 2 _18 - 65±0 03 for the SSM. The error bars represent 
the uncertainty due to different independent runs. Sim¬ 
ilarly, the values of c m are calculated from the limiting 
value of Pss,m as 8S(r) —*• 0 for r = 1. For the RRSM 
at p = 0.5 and for the SSM, the estimated values are 
c m = 0.60 ± 0.01 and 0.25 ± 0.01, respectively. The val¬ 
ues of d m are essentially the inverse of the averages of 
the scaled variable (8S(r)/r 2Hm ) for the respective mod¬ 
els. The values of d m are 2- 3 ' 66±0 02 and 2- 6 ' 93±002 for 


FIG. 10. (Colour online) (a) Plot of scaled toppling size dis¬ 
tributions P s ,mS Tsm of the two models RRSM at p = 0.5 and 
SSM against the scaled variable s/L Dam for different sys¬ 
tem sizes L = 256 (black line), 512 (red line), 1024 (green 
line), and 2048 (blue line), (b) Plot of Pss,mC 2Hm against 
8S(r)/r 2Hm for different values of r (see the legend) for 
p = 0.5 of the RRSM and the SSM for L = 2048. The same 
distributions given in (a) and (b) are plotted in (c) and (d), 
respectively, after normalizing by the respective non-universal 
metric factors. 

the RRSM at p = 0.5 and the SSM, respectively. In 
Fig. flOt'c') and flUIdb the scaled distribution functions 
are plotted against their scaled variables incorporating 
the metric factors of the respective models for both dis¬ 
tributions. It can be seen that a reasonable data collapse 
is observed in both cases. Hence, both the scaling func¬ 
tions f sm and g m for two models are essentially the same 
apart from the non-universal metric factors associated 
with them. Therefore, the critical state of the RRSM at 
p = 0.5 belongs to the so called stochastic or Manna uni¬ 
versality class, although the origins of such a universal¬ 
ity class in the two models are completely different. The 
stochasticity in the RRSM at p = 0.5 is due to the simul¬ 
taneous presence of two conflicting toppling rules (CTR 
and ATR) randomly at equal proportions, whereas the 
stochasticity in the SSM is externally imposed through 
the toppling rules. This is therefore independent confir¬ 
mation of the existence of the Manna universality class 
as it is observed in other models [2Cj in the context of 
APT on a diluted lattice. 


V. CONCLUSION 

A continuous crossover from RSM to SSM (Manna 
class) is observed in a random rotational sandpile model 
as the fraction p of lattice sites with the anti-clockwise 
toppling rule (and the rest of the lattice sites are with 
the clockwise toppling rule) varies from 0 (or 1) to 0.5. 
As p changes from 0 (or 1) to 0.5, the system passes 
through a series of non-universal stochastic models at 











each value of p. Finally at p = 0.5, at which there is 
maximum disorder in the toppling rule, the RRSM cor¬ 
responds to the Manna class. However, not only does 
the origin of stochasticity in the RRSM and SSM differ, 
but the macroscopic parameters identifying the critical 
steady states of these models differ significantly as well. 
A scaling theory for such a continuous crossover is de¬ 
veloped and verified numerically by estimating a set of 
critical exponents related to the avalanche properties as 
well as to that of the toppling surfaces. The values of 
the critical exponents satisfy all scaling relations among 
them for all values of p. This study then not only repre¬ 


sents a continuous crossover from the RSM to the SSM, 
but it also confirms the existence of the Manna class at 
the strong disorder limit of the RRSM. 
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